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Abstract 

The Lagrangian velocity structure functions in the inertial range of fully developed 
fluid turbulence are derived basing on the Navier-Stokes equations. For time r much 
smaller than the correlation time, the structure functions are shown to obey the 
scaling relations Kn{T) oc t^". The scaling exponents Cn are calculated analytically. 
The obtained values are in amazing agreement with the unique experimental results of 
the Bodenschatz group [I] . New notion - the Lagrangian position structure functions 
i?„(r) is introduced. All the Rn of the order n > 3 are shown to have a universal 
scaling. 

1 Introduction 

The spectrum of velocity pulsations in a turbulent flow can be naturally divided into three 
ranges: the large scales, the inertial range, and the viscous range. The largest scales L 
are comparable with the size of the turbulent volume. These are, roughly speaking, the 
scales of the largest vortices generated in the flow. The large-scale turbulence is a direct 
result of laminar flow instability for given boundary conditions. It is always non-uniform 
and anisotropic, the amplitudes of velocity pulsations are the largest at these scales. The 
most part of turbulent energy is contained in the large scales. On the other hand, viscous 
effects play a dominant role at scales smaller than the Kolmogorov scale r] = L/R^^^, where 
R is the Reynolds number [2]. As the viscosity decreases, the dissipative range becomes 
narrower. The intermediate range of scales ?7 <^ A <^ L is called the inertial range. In these 
scales, there is no energy dissipation; the energy flows in Fourrier space from large scales 
to the viscous range. Since viscosity is negligible in the inertial range, the Navier-Stokes 
equation can be reduced to the incompressible Euler limit [3]. 
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Hereafter we discuss the inertial range of locally homogeneous, locally isotropic, statis- 
tically stationary turbulent flow of incompressible liquid. The phenomenological theory of 
such turbulence constructed by Kolmogorov (K41 theory) predicts universal scaling laws 
for all statistical values. However, experiments have revealed departure from K41 caused 
by intermittency. 

In [4J we proposed a new model of turbulence explaining the intermittency and built on 
the Navier-Stokes Equation. Our idea was that the main role in statistics (in particular, 
in structure functions) belongs to the regions with very high vorticity. We showed that 
these regions stretch out taking the form of vortex filaments. We also showed that the 
growth of vorticity in the filaments is caused by large-scale pulsations. A considerable 
simplification of the problem is reached by treating the fully developed turbulence from 
Lagrangian viewpoint, which is natural from theoretical point of view [5]. We obtained the 
equation for the probability distribution function (PDF) of vorticity, and time dependence of 
its statistical moments. Note that recently the model equation interpreting the experimental 
data have been proposed [6]. This model is phenomenological and is not based on the 
Navier-Stokes Equation. 

This paper is a continuation of [1]. We apply the theory developed in [4] to obtain the 
Lagrangian statistical moments of velocity differences also called the Lagrangian structure 
functions 

k^{t) = {\v{t + t) - v{t)n , (1) 

where the values of velocity are taken along the trajectory of a liquid particle, and the angle 
brackets denote the average over the ensemble of trajectories. We find a relation connecting 
the Lagrangian structure functions of different orders. It may be used as an analogue to 
the extended self-similarity ansatz (established for the Euler case) [H [1] in the Lagrangian 
case. 

Recently there has been a significant experimental progress in measuring the Lagrangian 
statistical properties [H [8] . We compare the predictions of our theory with recent experi- 
mental results of the Bodenschatz group [1] and find a wonderful agreement. 

The paper is organized as follows. In §2, we recall briefly the method and some of the 
results obtained in [1] . In particular, we derive the time dependence of statistical moments 
of vorticity in a vortex filament. In §3, basing on these results, we calculate the scaling 
exponents of the Lagrangian velocity structure functions ([1]). We also derive a relation 
between the structure functions of different orders, which may be valid even in the cases 
when there is no scaling (for example, near to the boundaries of the inertial interval). In §4, 
the theoretical results are compared with the experimental data. In §5, we introduce a new 
notion: Lagrangian structure functions of position. We calculate their scaling exponents 
on the basis of our model. These results can also be checked by experiments. 
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2 Properties of turbulent structures in the inertial 
interval. 



We consider the Navier-Stokes equation for incompressible liquid. At the scales inside the 
inertial range the equation takes the form of the Euler equation: 

^ + (vV)v + ^ = 0, V-v = 0, (2) 

at p 

where v is the velocity of the flow and p is the pressure. The density p is taken to be unity 
below. The second equation expresses the incompressibility of the liquid. The restriction 
of scales by the inertial range presumes, in particular, that we consider only smooth initial 
conditions. The equation for pressure follows from ([2]): 

- Ap = ViVj ■ VjVi (3) 

Thus, the equations system ([2]) is complete. We note that ([2])can be rewritten as 

^=rot[v,a;], 

where w = rot v is vorticity. 

We now discuss the formation of stretched structures - vortex filaments - in the turbu- 
lent flow. They appear as a result of flow instability caused by large-scale pressure pulsations 
[1]. Random large-scale forces stretch a liquid drop forming a filament. The main role in 
the process belongs to the incompressibility of the liquid. First, the volume conservation 
results in transversal compression of the filament during its straining. Conservation of an- 
gular momentum then leads to acceleration of rotating and hence to growth of vorticity in 
the filament. Second, the sound speed in an incompressible liquid is infinite. This means 
an instant transport of large-scale pressure pulsations, which cause the stretching of the 
filaments. 

To illustrate these statements, we consider an axially symmetric flow. The hydrody- 
namic equations for the radial, azimuthal, and axial velocity components Vr, v^, and Vz 
are 

dVr ^ dVr dVr dp 

dt ^ dr ^ dz r dr 

^''+Vr^ + Vz^+'-^=0 (4) 



dt dr dz r 

dvz dvz dvz dp 

dt ^ dr ^ dz dz 
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We seek a solution of the system (jl]), ([5]) in the hnear form 

= u{t)r , Vr = a(t)r , Vz = h{t)z (6) 
The corresponding pressure following from (l6l) must be 

The equation then implies a relation between a and 6: 

2a + 6 = . (7) 

This means the fluid volume conservation. Indeed, we consider a cylindrical drop with a 
radius R{t) and a length Z{t) at an instant t. Then it follows from ([6]) that 

R = a{t)R, Z = b{t)Z, 

Integrating these equations, we make sure that the cylinder volume conserves: 

T:R{tfZ(t) = nRlZoexp [ (2a(ti) + b{ti)) dti = ttRIZq. 

Jo 

Combining ([6]) with (jlj), we obtain a system of ordinary differential equations: 

a + — u'^ = —Pi 

LU + 2auj = (8) 
b + b^ = -P2 

Differentiating the second equation of system ([8]) and substituting other equations, we 
obtain 

u = ~P2{t)uj (9) 

The function P2{t) has a meaning of pressure fall along the cylinder axis. Let us suppose 
that it is determined by large-scale pressure pulsations in turbulent flow. Hence, it is a 
complicated function of time; we assume that its time average is zero. Then the time 
intervals when P2{t) > and P2it) < are equally probable. However, at P2(t) > the 
function uj{t) oscillates, the oscillation amplitude changing weakly. On the contrary, at 
P2it) < 0, the function uj{t) grows exponentially. It is clear that the value of u grows 
on average. We note also that from ([8]) and ([7]) follows the proportionality of u to the 
cylinder length Z. Hence, such growth of u means stretching the cylinder. We also note 
that, despite the nonlinearity of the initial system (HI), (E]), the final equation ([9]) is linear. 
The nonlinearity acts only in the directions transversal to the cylinder axis. 

This example illustrates the behavior of a drop with large vorticity in a turbulent flow: 
it stretches out forming a filament, and the vorticity continues to increase. In |3] we have 
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analyzed the general equations ([2]) in terms of Lagrangian variables. We have shown that 
the development of filaments occurs in a way analogous to that in the considered example. 
We have found a linear equation analogous to ([9]) describing the growth of vorticity in the 
accompanying reference frame: 

Xn = —Pnk^k , Xi = Ui , Pnk = V„VfcP (10) 

Here Xj and P are the components of oj and the pressure in the Lagrangian reference 
frame. The matrix pnk is generated by large-scale pulsations. In order to obtain the 
statistical properties of the flow, it must be considered as a random quantity. Hence, ffTOl) 
is a stochastic equations set. 

Instead of three equations of the second order, we consider a set of six first-order equa- 
tions: 

Xi=yi yi = -PijXj (11) 
We introduce the joint probability density 

/(t,x,y) = (5(x-x(t))5(y-y(t))) . (12) 

Here x(t) and y{t) are the solutions of ( fTTI) at the given realization of pij and initial 
conditions; the average is taken over the ensemble of all possible realizations. The variables 
X and y are independent. 

Let the flow be locally homogeneous and isotropic. Then, assuming statistical indepen- 
dency of different components of the matrix pij, we obtain from (fTTI) the Fokker-Planck 
equation for / (see [1]): 



df ^ df 



dt dxk 



9y2 \ dy, 



(13) 



where t is normalized by the characteristic time of the probability density change that is of 
the order of the correlation time. 

For what follows we need an asymptotic time dependence of the moments < x^" >. To 
find it we consider the invariant moments of the even order 2P: 



n + m + k = P 



Integrating (fT3|) with appropriate weights, we obtain for these moments a closed set of 
linear differential equations 

^Mp(n, m, k) = 2nMp{n -l,m,k + l) + 2m{Ak + 2m + 2)Mp{n + l,m-l,k) 
+4m{m - l)Mp{n, m -2,k + 2) + kMp{n, m + l,k -1) + 2k{k - l)Mp{n + 2,m,k~2) 
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The number Q of equations in the set is equal to a number of combinations (n, m, k) such 
that n + m + k = P: 

Q = ^(P+l)(P + 2) 

Thus, the evolution of the 2P-order invariant moments is described by the set of Q linear 
differential equations. As t — > oo, the solutions increase exponentially, in particular (x^") oc 
exp{A2nt), A2n being the maximal root of the corresponding characteristic equation. Since 
< x^" > is positive, the solution must not oscillate, and the root must be real. To find the 
exponents, we calculate the determinants and solve the characteristic (algebraic) equations 
numerically. The first sixteen values A„ are 



A2 = 


= 2.52, 


A4 = 


= 6.12, 


A6 = 


= 10.43, 


A8 = 


= 15.25, 


Alo 


= 20.48, 


A12 


= 26, 


Al4 


= 32.03, 


A16 


= 38.25 


Ai8 


= 44.73, 


A20 


= 51.46, 


A22 


= 58.42 , 


A24 


= 65.58 


A26 


= 72.95, 


A28 


= 80.52, 


A30 


= 88.26 


A32 


= 96.16 



(14) 



To summarize, we cite [3] to list the main properties of the equation (fT3l) concerning the 
moments: 

1. All n-order moments of Xk and yj are connected by a set of first-order linear differential 
equations. Hence, it is possible to evaluate the moments of any order. 

2. The even moments grow exponentially. Independently of the initial conditions, the 
function / at large values of t depends only on the moduli x and y and on the cosine of the 
angle between the vectors /i = (x, y)/xy. 

3. The higher even moments grow faster than the lower ones. 

The properties 2 and 3 express the presence of intermittency. For example, for large 
values of t we have < x^" >^< x^ 

3 Lagrangian structure functions 

As shown in Section 2, regions of large vorticity in a turbulent flow take the form of filaments 
stretched along the vorticity direction. 

Then a typical trajectory in the flow demonstrates two types of behavior shown in 
Fig. 1. In the region I the vorticity is not very large, and the trajectory is smooth, the 
characteristic time of velocity change is of the order of the correlation time Tc- In the region 
II the vorticity is large, the test particle oscillates with frequency u. These second regions 
give the main contribution to correlation functions. Indeed, the second order Lagrangian 
structure function K2{t) is proportional to in "smooth" parts of trajectory, and we shall 
see below that in the "oscillating" parts it is proportional to r. 

So, we consider the regions with large u. We choose a cylindric frame with the origin 
in the local maximum of u, the axis z directed along u. Since there is fast rotation around 
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Fig.l. A typical trajectory of a test particle. 

the z axis near the origin, the flow does not depend on the angle 0. Then in the linear 
approximation we obtain the expression ([H]) which returns us to the example considered in 
Section 1: 

f0 = u;(t)r, t>r = a(t)r, Vz = b(t)z (15) 
Conservation laws for mass and for angular momentum along the z axis take the form (see 

dZD, m) 

2a + b = , 
uj = —2a{t)uj 

From ( fTSi) we find 

dr , . dz , , . dd) , , ,^ ^, 

- = a(t)r, ^=K^)^, = ^^^^ 

The parameters uj, a, h change slowly, with characteristic time t^. To the contrary, the 
oscillations are fast, uot^ ^ 1. Hence, with accuracy t/tc we have 

|5v| = |v(t + r) - v(t)| = 2w^sin Y = 2ra;sin^^ (17) 

To calculate the correlation functions of degree n we must raise an absolute value of 5w to 
the ra-th power and take an ensemble average with the probability distribution function. 
We first consider the pair correlation function 

K,{T) = {{v{t + r)-v{t)f) 

In accordance with the above, it has the form 

7^2 (r) oc / P(cj,t)tu^ sin^ (cur/2) rfcj (18) 
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In |1] we have analyzed the solution of (fT3|l and the integrated PDF P{x) = J f{x, y, fi)x'^dydfi 
(here fi is the cosine of the angle between x and y). We have found that P{uj) has an in- 
termediate stationary asymptote 

P{uj) = Cuj-^ (19) 

For small and for infinitely large values u, the PDF remains non-stationary. 

We see that the integral (fT8l) converges for the stationary PDF (fT9|) . Actually, 

r uj' sin^(W2) , rsin^ 

A2(t) oc / duj(XT I — dq 

Jo ^ Jo Q 

(The input of Sv"^ is small - proportional to z'^). 

So, for T <^ Tc the second order structure function along the Lagrangian trajectory is 

K2{t) cx r (20) 

This result adjusts with the Kolmogorov's theory K41 [2], [2]. 
For higher-order structure functions Kn{T) we have 

POO 

Js:„(r) oc / P(cj,t)cj"|sin(c<;r/2)|"c/a; (21) 
Jo 

The integral with the stationary PDF diverges. Hence, the higher structure functions are 
determined by the non-stationary part of the PDF. To calculate i^„(r), we use the moments 

P{uj,t)uj^''duj oc e^2"* 

found in Section 2. We expand the sine in fl2T|) into the series and integrate each term. For 
time t satisfying 

ln(r-2^ 



1 << t << ^ — , r^O 



the second item is much less than the first one. We restrict ourselves by the first term as 
T 0. Then 

ir„(r) oc r"e^2"* (22) 

Excluding t, we obtain 

'i^™(r)V^^''" /C^nl^^i^n(r)^'/^^" 



r" 



(23) 



Here Cn\m are constants not depending on r. This relation connects the structure functions 
of different orders. Choosing m = 2 and taking ( !20l) into account, we find 



K„, oc r' 



Cn 
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where 

Cn = n-^ (24) 
Using (fT^ we get first twelve scaling exponents (n- 

Ci = 0.59, C2 = l, C3 = 1.3, (4 = 1.51, 

(5 = 1.65, C6 = 1.75, C7=1.77, C8 = 1.T5, 

Cg = 1.69 , Cio = 1.59 , Cii = 1.45 , C12 = 1.28 



Actually, the stationary structure function ( ]20l) is not necessary to determine the scaling 
exponents. The reason is that the parameters grow fast as a function of n. For example, 
we take m = 16. Using (|23|1 with A32 = 96.16, we then obtain 

Knir) oc r"-^^^2"/^32 K^i"^^'' 

For small n the ratio h.2n/ is small enough, and we can put i^j^^"^^^^ ^ 1. For n = 
2 we then find K2{t) oc t^'^^ ~ r. The result is very close to that obtained with the 
stationary PDF. We note that this relation is derived now from ( l23i) only. This demonstrates 
the consistency between the stationary and non-stationary solutions of the equation (fT3|) . 
Taking other values of n, we get in the same way 

C( = 0.58, C2 = 0.98, C3 = 1.26 

C: = 1.46, C5 = 1.59, C6 = 1.67 

We recall that in the derivation of ( l23i) we neglected the rest of the sine expansion in ( l2Til 
as r — J. 0. To justify this approximation, we calculate ( l2Ti) with complete sine expansion, 
using the moments of u. We then express e* through Kiq using (1221) . For example, for K2 
we have 

K2 = -T^e^^' - — r^e^«* + ^r^e^^* + ... (25) 
4 48 1440 ^ ' 



_ 1 ^0.98 jv-O.Oe 1 2.26 tv-O.II , ^ -3.46 jv'0.16 , 

-4^ 48^ + 1440^ 

We see that as r — > the series converges rapidly, and the rest is small with respect to the 
first term. The same is correct for the structure functions of higher orders. 
We note also that the relation fl23|) can be rewritten as 



7>'A2n/A2m . 



2n 



Cn\ra = 77 ^ ^"^ (26) 

According to our theory, these combinations of structure functions must be constant, i.e. 
they must not depend on r. 
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4 Comparison of the theory with the experiment. 



In this Section we compare the predictions of our theory with the results of recent exper- 
iments of the Bodenschatz group [H],II]- In these experiments turbulence was generated 
by two counter-rotating disks. Different regimes of the facility allowed to make the mea- 
surements for different local Reynolds numbers Rx. The Lagrangian test particles used 
to trace the flow had the size comparable or smaller than the Kolmogorov scale rj for all 
tested Reynolds numbers. The motion of particles was tracked in a subvolume of about 
10~^ of the whole turbulent volume, in the center of the tank where the effects of the mean 
velocity were negligible. Thus, this region seems to us optimal for investigation of locally 
homogenous turbulence. 

The measurements performed in [8] showed directly the presence of stretched structures 
- vortex filaments. It was shown that near to these structures the absolute values of particle 
accelerations were maximal. The probability distribution of these accelerations was found 
to be significantly non-gaussian in the region of their large values. Thus, general picture of 
vortex structures in turbulence corresponds to that in our theory [1]. 

Table 1. Values of the scaling exponents normalized by ^2- The experimental results for 
different local Reynolds numbers are cited from [Ij 



R 


200 


690 


815 


Theory 


C1/C2 


0.59 ±0.02 


0.58 ±0.05 


0.58 ±0.12 


0.59 


C3/C2 


1.24 ±0.03 


1.28 ±0.14 


1.28 ±0.30 


1.3 


C4/C2 


1.35 ±0.04 


1.47 ±0.18 


1.47 ±0.38 


1.51 


C5/C2 


1.39 ±0.07 


1.61 ±0.21 


1.59 ±0.46 


1.65 


C6/C2 


1.40 ±0.08 


1.73 ±0.25 


1.66 ±0.53 


1.75 


C7/C2 


1.39 ±0.09 


1.83 ±0.28 


1.67 ±0.60 


1.77 


C8/C2 


1.40 ±0.10 


1.92 ±0.32 


1.65 ±0.66 


1.75 


C9/C2 


1.42 ±0.11 


1.97 ±0.35 


1.61 ±0.73 


1.69 


C10/C2 


1.46 ±0.12 


1.98 ±0.38 


1.57 ±0.80 


1.59 



We make the quantitative comparison of the Lagrangian structure functions obtained in 
our theory with the results of the experiment [T] . Table 1 and Fig. 2 show the Lagrangian 
scaling exponents (n normalized by ^2, up to the 10th order. The first three columns of 
Table 1 represent the values measured in [1] with different local Reynolds numbers Rx. The 
last column contains our theory prediction (^M . The theory is in excellent agreement with 
the experimental data. Moreover, the agreement is better for larger values of Rx- This 
corresponds to the fact that the theory is constructed in the limit R — > 00. Also, both in 
the theory and in the experiment with the most value of Rx, the scaling exponents grow 
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up to the 7-th order, and then decrease. We note that the theory predicts also the scahng 
exponents of higher orders that are not measured yet. 




c I I I I I I I I I I I 

1 23456789 10 

n 

Fig. 2 Scaling exponents Cn normalized by C2 as a function of order. Small black squares denote 
the prediction of our theory, other symbols show the experimental results [1] for different local 
Reynolds numbers: Rx = 200 (red squares), Rx = 690 (green circles), and Rx = 815 (blue 
triangles). The solid line corresponds to the Kolmogorov theory K41. Strong departure from the 
K41 prediction of both the experiment [1] and our theory is the manifestation of intermittency. 

Table 1 and Fig. 2 also demonstrate the growth of dispersion of the measured scaling 
exponents as a function of order and of Rx- This is in accord with the presented theory, 
which is based on the notion that filaments (narrow regions with very high vorticity) make 
the most contribution to structure functions. The dispersion behavior is then the result of 
intermittency. The higher is the order, the more important become very high peaks, the 
more seldom they occur. Hence, for a given sample the dispersion increases with the moment 
order. On the other hand, the real height of peaks is limited by viscosity]^ As viscosity 
decreases, the size of a sample needed to determine the statistical moments increases. This 
property of statistical systems is demonstrated, in particular, in [9]. 

^Presumably, the viscosity constrains the value of m in (|23p . We choose m = 16 as an upper limit 
because the ratio of the large-scale characteristic time Tc to the viscous Kolmogorov time is about 100. 
Since the time of change of Kiq is Agg^ O.OItc, it lies near the boundary between the inertial and the 
viscous ranges. 
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In addition to analysis of the scaling exponents, we also compare (l23l) with direct ex- 
perimental data [1] for dependence of the high-order Lagrangian structure functions of r. 
For that we choose the form ( !26l) of the equation (!23l) to check if the combinations C„|m of 
any two structure functions are constant. The results for Cio|9 and Cio|8 are presented in 
Fig. 3. We see that even for the highest of the measured orders, the values C are constant 
up to the accuracy of the available experimental graph. 



1 E1 4 — I 1 1 — I — I — I — r-| 1 1 1 1 — I — I — I I I 1 1 1 1 — I — I — I — r 




1 10 100 



t/t 

n 

Fig. 3. Lagrangian structure functions of the orders 8,9, 10 (black curves from the bottom to the 
top) normalized by the K41 predictions {R\ = 690): cited from [1]. Red and green curves denote 
the functions Cio|9 and Cio|8) respectively. 

We stress that no adjustable parameters were used in calculating both the scaling expo- 
nents (Table 1, Fig.2) and the combination of the Lagrangian structure functions (Fig. 3). 
Thus, the coincidence between the theoretical predictions and the experimental results is 
wonderful. One can suppose that a significant role belongs to a happy choice of parameters 
of the experiment: very homogeneous flow in the subvolume where the measurements were 
performed, and the size of test particles coinciding with the Kolmogorov scale. Due to 
this the high frequency noise is damped in the measurements, and the experimental data 
describe the inertial range only. Exactly this range is studied by the theory. 
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5 Universal scaling for Lagrangian structure functions 
of position. 

The Lagrangian approach gives a possibihty to study not only the statistical properties of 
velocities, but also that of positions of the particles. This allows to introduce a new notion 
in the turbulence theory. Define the Lagrangian structure functions of positions of two 
particles as 

i?„(r) = (|r(t + r)-r(t)r) 

As in the case of velocity structure functions, the main contribution in Rn comes from the 
regions with large vorticity. Integrating (fT6|) and calculating \6r\ by analogy with (fT7|) . we 
obtain 

K(r)oc J P(c^)sin"(c^r/2)rfcuoc j ^^^"(^^^/^) 

If n > 3 the integral with the stationary PDF f|T9l) converges. Hence, the high-order 
structure functions of positions are determined by the stationary part of the PDF. We note 
that in the asymptotics all the scaling exponents are equal: 

/■ sin"(cjr/2) , o 
Rn(r) oc / , du} (xr\ n > 3 



In other words, there is a universal scaling for all n > 3. This universal scaling, if it 
would be observed in experiments, can be an evidence for the determining influence of the 
filaments on the behavior of the structure functions in the inertial range. 



6 Conclusion 

To conclude, we list the main results of the paper. 

1. The Lagrangian velocity structure functions of high orders are derived in the inertial 
interval of fully developed liquid turbulence. 

2. The comparison of the theory predictions with experimental results shows an excellent 
agreement. 

3. New experimentally measurable values - the Lagrangian position structure functions 
are introduced and calculated. Their scaling exponents are shown to be independent 
on the order of the function for any n > 3. 

We emphasize that no adjustable parameters were used for the scaling exponents. Also 
the power-law dependence of the Lagrangian structure functions was not suggested in the 
model. All the obtained relations are the consequences of the equation f|T3|) derived directly 
from the Navier-Stokes equation in the inertial interval. 
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